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RIGIDITY SPECTRUM OF FORBUSH DECREASE 
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ABSTRACT 

Using data from neutron monitors and muon telescopes at 
surface and underground stations (Nagoya , 0 m.w.e.; Misato, 

34 m.w.e.; Sakashita, 80 m.w.e.), we obtain the average 
rigidity spectrum of Forbush decreases (Fds) during the 
period of 1978-1982. Thirty eight Fd-events are classified 
into two groups "Hard Fd" and "Soft Fd" according to size 
of Fd at Sakashita station. It is found that a spectral 
form of fractional-power type (P“ Y 1 (P-fP c )”Y 2 ) is more suit- 
able for the present purpose than that of power-exponential 
type or of power type with an upper limiting rigidity. 

The best fitted spectrum of fractional-power type is ex- 
pressed by y ,=0.37 , y 2 =0.89 and P C =10GV for Hard Fd and 
y 1 =0.77, y 2 “1*02 and P C =14GV for Soft Fd. 

1. Introduction. Rigidity spectrum of Forbush decrease (Fd) has been 

studied by many authors, using neutron monitor data(e.g. Lockwood, 1971; 
Iucci et al., 1979; Lumme et al. , 1983; Fenton et al., 1983), or by 
adding a small amount of underground muon data to those in the above (e.g. 
Thambyahpillai et al. , 1965; Andreis-Sandrucci et al., 1968; Mishima et 
al. , 1973; Suda et al. , 1977; Krymsky et al., 1979; Wada et al. , 1979; 
Sakakibara et al., 1979; Wada and Suda, 1980). These studies have been 
made by assuming apriori a definite form of the spectrum and have tried 
to determine parameters involved in the spectral form. In the present 
paper, we make a comparative study of the spectral forms of Fd, using 
plenty of data in wide range of primary cosmic ray rigidity. For the 
muon telescope data, we use the response function given by Murakami 

et al. (1979,1981) . For the neutron monitor data, we choose the re- 
sponse function given by Nagashima (1980) for the maximum solar activity 
out of three functions (Lockwood and Webber, 1967; Aleksanyan et al.,1981; 
Mori and Nagashima, 1984) , based on a comparative study of response func- 
tions (S • Sakakibara et al.,1984; S.Mori and K. Nagashima, 1984) . 

2. Rigidity spectrum of Forbush decrease. Data used for the present 

analysis are those of 8 neutron monitors (median primary rigidity P m 21 ^ 
35GV) , 17 muon telescopes at Nagoya(P m =60 ^ 119GV) , 9 underground muon 
telescopes at Misato (P m = 14 5 ^ 209GV) and 3 underground muon telescopes 
at Sakashita(P m =330 ^ 567GV) during the period of Feb. 1978 - Feb. 1983* 
Fourty five Fds greater than 2% in depression for neutron intensity at 
Alert station were selected, but because of missing data of muon compo- 
nents, 38 Fds were finally adopted. 

Fig.l shows the rigidity dependence of the development of Fd on daily 
basis, obtained by superposing those 38 Fds at the epoch which 0-day in- 
cludes the occurence time of SC. As can be seen in this figure, the 
depression at Sakashita is about 0.15% on the average and statiscally 
significant, but individual Fd at the station is not always statiscally 
significant and, in some cases, shows some increase. Taking this fact 
into account, we classify the Fds into two groups "Hard" (23 events) and 
"Soft" (15 events) according to size of Fd of Sakashita vertical component. 
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The magnitude(AI) of each Fd is 
defined by the difference between 
daily mean intensities on +1 day 
and -3 day(cf. Fig.l). The bound- 
ary of two groups is set at the 
magnitude(AI=0.05%) of Sakashita 
vertical component, which corre- 
sponds to the statistical error 
for the component. 

We made a comparative study 
of the rigidity spectrum of Fd, 
using the following three types 
of the spectrum. 


PW : Power type 


AI<*{ 



for P < P u 
for P > P u 


PE : Power-Exponential type 
AI P Y exp(-P/P 0 ) 


Fig.l Averaged daily development of 
38 Forbush decreases, superposed at 
the epoch day including SC, and 
that of the corresponding solar 
wind velocities. 


FP : Fractional-Power type 
AI « P _Yl (P+Pc)~ Y * 

( 1 ) 


Corresponding to the parameters involved in Eq.(l), the expected values 
Alexp* s f° r m uon components and neutron monitors are calculated, respec- 
tively, from the differential coupling coefficients by Fujimoto et al. 
(1984) and from those by Mori and Nagashima(1984) . By the least squares 
method, we found that FP-type spectrum in Eq.(l) can produce the expected 
values AI e xp fs most compatible with the observations. To demonstrate 
the coincidence betwee'n the observation and AI eX p f s, we use a diagram 
whose ordinate represents the observed magnitude of Fd and whose abscissa 
expresses the corresponding magnitude Al eX p calculated by using a rigidity 
spectrum of Fd. We call this the correlation diagram of Fd. If AI e xp 
is calculated by using the spectrum determined from the observed data by 
the least squares method, we call the diagram the best fitted correlation 
diagram of Fd. In this case, all points on the graph must lie on a. line 
inclined by 45° from the abscissa. Such diagrams of Hard and Soft Fd’s 
are shown in Fig. 2, for the FP-type spectrum. It is noted that, in case 
of Soft Fd, the least squares method was applied to those data excluding 
those of Sakashita station. 

The rigidity dependence of Fd’s can be expressed more plainly if we 
use an effective primary rigidity(P]7) for Fd, which is defined as the 
median rigidity of a function obtained by multiplying the response func- 
tion by the best fitted rigidity spectrum of the Fd. Fig. 3 shows the 
Pg-dependence of Hard and Soft Fd T s. The average magnitudes <AI>’s can 
line up systematically on this graph. 

In Table I, we show the best fitted parameters of three types of the 



Average Magnitude of Fd <AI> 
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Table I. The best fitted parameters for three types of the 
rigidity spectrum. 


\Spectral 

^\typ e 

Group^^^ 

rr p _Ti (p+p c ) Y2 

PE 

P ^exp (-P/Po ) 

ru ? Y , (P « P„) 

1 Y i Y 2 P C (GV) 

Y Po (GV) 

P U (GV) 

Hard Fd 

0.37 0.89 10 

0.90 471 

0.98 400 

Soft Fd j 

0.77 1.02 14 

1.17 145 

1.26 100 



Fig. 2 The best fitted correlation 
diagram of Hard and Soft Fds. 
AI exp is calculated by using 
the fractional-power type spec- 
trum with parameter shown in 
the figure, and in case of neu- 
tron monitors the NAGASHIMA- 
response function is selected 
for the calculation. 


I i m T] 1 1 1 — r 1 T'T rr- 

“T"' 1 1" "1 

I 8 Neutron st. 

- 

- Nagoya 


b 1 i Mi sato 4 

1- 1 

Sakashita -j 



Effective Primary Rigidity (P e ) GV 


Fig. 3 Pg -dependence of average mag- 
nitude <AI> of Hard and Soft Fds. 
The effective primary rigidity Pg 
for Fd defined in the text is cal- 
culated under the same conditions 
as i n Fig. 2. Errors are derived 
from the dispersion of individual 
Fds. 
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3. Conclusion . It is concluded that the FP-type spectrum is more 

suitable than PE-type or PW-type to express the observed rigidity depen- 
dence of Forbush decrease. However, in case of the Soft Fd, one cannot 
select the best spectral form out of the three types, by reason of the 
insufficient accuracy of the observed data in the higher rigidity range. 
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